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Thermoelectric properties of the system La2Ni04+<5 have been recently discussed [Phys. Rev. B 
86, 165114 (2012)] via ab initio calculations. An optimum hole-doping value was obtained with 
reasonable thermopower and thermoelectric figure of merit being calculated. Here, a large increase 
in the thermoelectric performance through lattice strain and the corresponding atomic relaxations is 
predicted. This increase would be experimentally attainable via growth in thin films of the material 
on top of different substrates. A small tensile strain would produce large thermoelectric figures of 
merit at high temperatures, zT ^ 1 in the range of oxygen excess 5 ~ 0.05 - 0.10 and in-plane lattice 
parameter in the range 3.95 - 4.05 A. In that relatively wide range of parameters, thermopower values 
close to 200 fiV/K are obtained. The best performance of this compound is expected to occur in 
the high temperature limit. 



I. INTRODUCTION 



The thermoelectric (TE) performance of a mate- 
rial is usually quantified by means of the so-called di- 
mensionless figure of merit zT= <jTS 2 /^, where S is 
the thermopower, a the electrical conductivity and 
k the thermal conductivity that can be expressed 
as the sum of electronic and lattice contributions, 
zT > 1 is required for applications^ 
and in order to maximize zT, a high S, high a 
and low k are required. All the electronic mag- 
nitudes (electrical conductivity, thermopower and 
the electronic component of the thermal conductiv- 
ity) are inter-related and their simultaneous opti- 
mization is somehow conflicting. Low-gap semicon- 
ductors are the best meeting the compromises con- 
sidering only the electronic properties. If one re- 
duces (increases) the band gap of a semiconductor, 
the conductivity increases (decreases) and the See- 
beck coefficient goes down (up). Thus, there is a 
compromise situation at intermediate doping levels 
where TE efficiency is enhanced. To reduce the lat- 
tice thermal conductivity, nanostructuring is a com- 
mon strategy. 3 The situation gets more complicated 
when instead of standard parabolic-band semicon- 
ductors, strongly correlated electron systems with a 
low band gap are considered^ These can in princi- 
ple be equally effective from an electronic structure 
point of view, but they provide additional ingredi- 
ents that allow to tune the band structure in order 



to optimize separately the different magnitudes in- 
volved in the TE response: having localized elec- 
trons allows for the tuning of band splittings, band 
widths, correlation and ordering effects, etc. More- 
over, band engineering is possible in strongly cor- 
related electron systems in order to tune the elec- 
tronic structure for an enhanced TE response, i.e. 
they are more flexible than standard semiconduc- 
tors for a tunability of their band structure through 
magnitudes other than doping or nanostructuring 
(commonly used for sp semiconductors). Parame- 
ters such as strain, ordering phenomena, electron 
localization mechanisms, pressure effects, etc. often 
play a big role in their electronic structured In par- 
ticular, strain engineering has recently been shown 
as a direct method to control the TE properties in 
topological insulator materials Bi2Se3 and Bi2Te3. 6 

Among oxides and other 3d electron systems, mis- 
fit layered cobaltates have drawn significant atten- 
tion as TE materials 1,7 10 because of the large ther- 
mopower observed together with metallic conduc- 
tivity. Their interesting TE properties stem from 
the existence of two electronic systems within those 
cobaltates. 11 Due to the trigonally distorted octa- 
hedral environment, the Co t2 P bands are split be- 
tween a highly mobile wide band that drives metal- 
lic conduction (of e* character) and a narrower band 
(with aip parentage) that provides the large Seebeck 
coefficient found in that family of oxides. 12 Similar 
electronic structure and promising TE properties via 



hole-doping have been predicted in La2Ni04+5 by 
means of ab initio calculations^ The existence of lo- 
calized and itinerant electrons in that system is also 
at the heart of a reasonable TE response. In that 
paper, the electronic structure under doping was an- 
alyzed in order to understand which bands are influ- 
enced by doping and what amount of oxygen excess 
is needed to enlarge the TE efficiency. Experimen- 
tally, the relevant TE parameters observed typically 
at room temperature are S ~100 /iV/K, 14 p ~ 5 
mll.c m 15 ' 16 and n ~ 7.5 W/m.K^ 7 - which add up to 
a zT at room temperature of about 0.01. This is 
comparable to other promising compounds such as 
CrNilSr-ffl 

Ref. [H suggested that the appropriate oxygen 
content 5 together with a thin film growth could 
significantly enhance zT. Growing thin films in- 
troduces the additional ingredient of strain due to 
lattice size mismatch when the oxide is grown on 
different substrates. In this paper we will analyze 
how one can refine even further the TE properties 
of La2Ni04+# when both oxygen-excess and strain 
are introduced. We will perform a band engineering 
study of how these strain effects can enhance the TE 
performance in a 3d electron system that has been 
predicted in the past to yield an interesting TE re- 
sponse. The goal of this work would be designing 
the proper growth conditions for an improved TE 
efficiency. 

The paper is organized as follows: Section HI1 will 
present the calculation details, Section IIIII will de- 
scribe the electronic structure and its strain depen- 
dence, analyze the TE properties calculated for var- 
ious oxygen contents, and finally we will summarize 
the main conclusions of the work in Section IIVI 



II. COMPUTATIONAL PROCEDURES 

Our electronic structure calculations were per- 
formed within density functional theor y 21 ! 22 using 
the all-electron, full potential code wien2k 23 based 
on the augmented plane wave plus local orbitals 
(APW+lo) basis set^ For the exchange-correlation 
functional in the structural relaxations we have used 
the Perdew-Burke-Ernzerhof version of the general- 
ized gradient approximation^ (GGA). To study the 
effects of strain, we have analyzed different in-plane 
a,b lattice parameters, and for each of them the c pa- 
rameter was calculated. Consequences will be drawn 



for strains in thin films but the actual calculations 
are done in bulk-like unit cells with tetragonal sym- 
metry for various c/a ratios. 

For the calculations of the transport proper- 
ties, we have used the recently developed Tran- 
Blaha modified Becke- Johnson (TB-mBJGGA) 
potential. 26 This has been shown to provide an accu- 
rate account of the electronic structure of correlated 
compound s 20 ! 27 using a parameter-free description, 
without invoking for a particular value of U selected 
by hand as typically done in the LDA+U method, 
and also yielding accurate band gaps for most semi- 
conductors. 

The transport properties were calculated using a 
semiclassical solution based on Boltzmann's trans- 
port theory within the constant scattering time ap- 
proximation by means of the BoltzTrap code, 
that uses the energy eigenvalues calculated by the 
wien2k code. 28 We refer the reader to Ref. [51 for 
details on how the different transport coefficients 
are obtained. In this case, denser k-meshes are re- 
quired, in our case up to 40 x 40 x 15 to reach 
convergence. The constant scattering time approxi- 
mation assumes the relaxation time r(e) as energy- 
independent. This results in expressions of both 
the thermopower and the TE figure of merit with 
no dependence on r (they can be directly obtained 
from the band structure without any assumed pa- 
rameters). This approximation has been used suc- 
cesfully to describe several potentially useful TE 
materials *22r— 

The simulation of different doping levels was per- 
formed in two different ways: i) by a shift in the 
chemical potential for the undoped compound, cal- 
culating the hole concentration at each temperature. 
We have seen in previous worksi^ that some of the 
essential features and trends can be obtained by this 
method, specially at extremely low doping, 27 ii) us- 
ing the virtual crystal approximation (VCA), that 
provides a more accurate description of the doped 
material and yields substantial differences in quanti- 
tative values of the TE properties, as we will discuss 
in detail below. Both calculations treat the dopants 
in an average way, whereas the actual oxygen-excess 
atoms will reside in some particular positions in the 
lattice and can lead to polaronic effects, dopant or- 
derings, etc. that will be partly missed by our cal- 
culations. In the low-doping regime our calculations 
will be more reliable. 

All calculations were fully converged with respect 
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to all the parameters used. In particular, we used 
R m tK max = 7.0, a k-mesh of 10 x 10 x 4, and muffin- 
tin radii of 2.35 a.u. for La, 1.97 a.u. for Ni and 1.75 
a.u. for O. 



III. RESULTS 

A. Revisiting the electronic structure and 
magnetic properties of unstrained La2NiC>4 

Let us recall some basic ideas about the electronic 
structure of La2Ni04 + s. This is a compound con- 
taining Ni 2+ :d 8 cations in a largely elongated oc- 
tahedral environment. The crystal field splittings in 
such geometry lead to a breakup of the e^ degeneracy 
with the d x 2_ y 2 band higher in energy than the d z 2 
(as sketched in Fig. QJi). A significant Hund's rule 
coupling stabilizes the high-spin state (HSS) t2g 6 e g ° 
(S= 1). In previous works j 13 i 34 it was seen that 
hole-doping via oxygen exccess introduces holes in 
the Ni d x 2_ y 2 band, which is somewhat localized. 
This d x 2_ y 2 band below the Fermi level can be sub- 
stantially modified in a process one can call "band 
engineering" in order to improve the TE response 
of the compound. To that end, one needs e.g. a 
larger thermopower (proportional to the derivative 
of the density of states (DOS), and hence increas- 
ing with a reduced bandwidth) without compromis- 
ing the electrical conductivity. One possibility for 
that is to bring the d z 2 band (whose bandwidth is 
smaller due to the comparably smaller off-plane hop- 
ping) closer to the Fermi level, and make use of its 
reduced band width to enlarge the thermopower. Of 
course, this is not exactly so because, as mentioned 
above, the electronic-only magnitudes in semicon- 
ductors are all closely related, and in principle a re- 
duced band width might lead to an increased band 
gap and a reduction of the conductivity. Yet an op- 
timum performance can be sought, as we will see 
below. In any case, to bring the d z 2 higher in en- 
ergy (closer to the Fermi level for the majority spin 
channel), this band needs to be destabilized with re- 
spect to the d xy (which itself is split from the lower- 
lying d xz /d yz doublet) and d x 2_ y 2 bands (see Fig. 
[1]). Thus, enlarging the a,b plane could produce the 
required effect. This would, however, significantly 
reduce the in-plane conductivity, which is the largest 
component in this layered compound, with the cor- 
responding fall-out in figure of merit (let us recall, 



zT=aS 2 T/ k). If, however, a,b lattice parameters 
are reduced, a can be increased, but the d x 2_ y 2 band 
will be moved up in energy so that it will remain the 
only accessible band via hole-doping (a sketch of how 
all these bands shift with strain can be seen in Fig. 
[I]). In that case, the band will become broadened 
(lower peaks for the DOS, smaller thermopower) by 
the increased in-plane hopping caused by the smaller 
lattice parameter. Hence, there will be optimum 
values (a compromise) for both strain and doping 
(trying to tune the oxygen excess to reach the chemi- 
cal potential where thermopower is maximized) that 
one could use for designing thin films with optimal 
conditions so that an enhanced TE response can be 
obtained. Changes in the in-plane lattice parame- 
ters are induced by growing epitaxial thin films on 
different substrates. As explained above, a thin film 
geometry is beneficial for the overall TE response 
since k will be drastically reduced and a has been 
observed to increase due to its large two-dimensional 
anisotropy. 15 

Concerning the magnetic order, it has been shown 
that La 2 Ni0 4 (with Ni 2+ :d 8 S=l, HSS) is an anti- 
ferromagnet with an in-plane ordering such that an 
antiferromagnetic (AF) interaction between nearest 
neighbor Ni atoms is stabilized through an e^ 2 -0- 
e g 2 superexchange interaction. The commensura- 
bility of the ordering has been, however, put into 
question. 35 Careful studies with respect to the oxy- 
gen content show how the Neel temperature is sig- 
nificant even at values of 5 ~ 0.1^ and could vanish 
at S ~ 0.14J^ We will focus our calculations in that 
doping level region (5 < 0.15), where magnetic or- 
der exists. All our calculations assume an in-plane 
checkerboard AF ordered phase. 



B. Changes in the electronic structure caused 
by strain 

We have performed calculations at various in- 
plane lattice parameters, that try to simulate both 
compressive and tensile strains caused by growth of 
La2Ni04+£ thin films on top of different substrates. 
Calculations are carried out in bulk-like unit cells 
but help us draw conclusions about the changes in 
the electronic structure undergone by the compound 
if grown epitaxially on top of different substrates, 
with different in-plane strains and c/a ratios. 

We have optimized the out-of-plane c lattice pa- 
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FIG. 1: (Color online.) Sketch of the crystal field levels of a Ni 2+ :d 8 cation in a HS state to observe the effects of: 
a) a distortion below octahedral symmetry in an elongated octahedron, and the effects of both compressive (b) and 
tensile (c) strain. The gap is expected to be reduced (increased) for the compressive (tensile) strain case. 



rameter of the structure for each value of the in- 
plane a, b lattice parameters chosen. We have per- 
formed GGA calculations for several values of the 
in-plane lattice parameters from 3.7 to 4.1 A. The 
bulk value of the a lattice parameter in the tetrago- 
nal structure^ 7 of La2Ni04 is 3.89 A. For each a, c 
combination, the atomic positions were relaxed and 
the minimum-energy c-value was obtained for each 
in-plane area. In Table U we can see the results of 
the lattice parameter optimizations. Seeing the bulk 
value of the lattice parameter for La2Ni04, it is clear 
that when grown on a SrTiOs substrate, it would be 
roughly unstrained. Films have been grown in the 
past 15 on NdGa03 (with slightly lower lattice pa- 
rameter, somewhat similar to LaAlOs). Other sub- 
strates can be found with different lattice parame- 
ters, we have tried to study the evolution of the prop- 
erties with strain giving an interval of a,b in-plane 
lattice parameters that covers most of the typically 
available substrates. Probably the extreme cases we 
have considered will not be stable in the laboratory, 
but still help us drawing conclusions about the evo- 
lution of the transport properties with strain. 

TABLE I: In-plane vs. out-of-plane optimized lattice 
parameters for various strains. A compression in the 
xy-plane leads to an elongation in c to try to keep the 
volume constant. 



a (A) 


c(A) 


3.74 


9.39 


3.85 


8.98 


3.93 


8.85 


4.00 


8.72 


4.11 


8.51 



We know from previous studies and the discus- 
sions above that the only band available in un- 
strained La2Ni04 + 5 to be populated with holes via 
oxygen excess is the majority-spin d x 2_ y 2 band. 
Compressive strain (reducing the in-plane lattice pa- 
rameter) will move this d x 2_ y 2 band further up in 
energy but at the same time it will increase its band- 
width (together, an increase in conductivity and de- 
crease in the derivative of the DOS, and hence a 
reduction of the thermopower, is expected). If, on 
the other hand tensile strain is applied, this band 
moves further down in energy (it gets stabilized by 
the larger in-plane lattice parameter) closer to the 
energy window where the majority d z 2 band and the 
minority t2 g bands reside. Thus, the effects of ap- 
plying a tensile strain are slightly more complicated 
because more bands are involved. Moreover, tensile 
strain would reduce the in-plane conductivity (due 
to the larger in-plane lattice parameter) but could 
enlarge the Seebeck coefficient if a large DOS is re- 
tained close to the Fermi level (accessible via hole- 
doping). 

Figure [2] shows the band structure for the two 
types of strain studied: compressive and tensile. We 
see that compressive strain (on the left) reduces sig- 
nificantly the band gap, driving the system towards 
a bigger in-plane conductivity caused by the reduced 
Ni-Ni in-plane distance and corresponding increase 
in the Ni d x 2_ y 2 bandwidth (larger in-plane hopping 
mediated by those orbitals via a large a-bond with 
oxygens). The d x 2_ y 2 band is significantly less wide 
in the right panel (tensile strain) due to the reduced 
hopping caused by the enlarged in-plane lattice pa- 
rameter. Moreover, tensile strain leads to a stabi- 
lization of the d x 2_ y 2 band, that places it closer to 
all the other occupied bands (we see the band very 
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FIG. 2: (Color online.) Band structure (calculated within GGA) for two types of strain: compressive on the left 
(a= 3.74 A) and tensile on the right (a= 4.11 A). Notice the displacement of the d x 2_ y 2 band (easy to identify with 
maxima at X and U, being 1 eV wide, and occuring just below the Fermi level specially in the left panel) to a lower 
position in the tensile strain case. A significant reduction of the gap occurs for the compressive strain case and a 
somewhat smaller band width of the d x 2_ y 2 band is obtained in the tensile strain case. 



distinctly below the Fermi level in the compressive 
strain case, where the opposite is true, but some- 
what mixed with many other bands for the tensile 
one). 

There are various factors at play here that can 
enhance or reduce the different TE properties: con- 
ductivity and thermopower. Thus, calculations need 
to be performed to account for them all (doping, 
strain, changes in bandwidth, introduction of more 
bands near the Fermi level, etc.) properly. All in all, 
we have various conflicting parameters to tune and 
calculations are required to see the balance between 
them in order to obtain a compromise that might 
enhance the TE response. 



C. Influence of strain on the thermoelectric 
properties 

We can quantify all the important magnitudes for 
the TE response of a material in terms of the dimen- 
sionless TE figure of merit zT= ctTS 2 /k. As men- 
tioned above, zT needs to be larger than 1 for appli- 
cations. It is often quoted that a thermopower larger 
than 200 /xV/K is needed for a high-performance 
TE response. We can, from first principles calcula- 
tions through the electronic structure of the mate- 
rial, give an estimate of the electronic part of the 
figure of merit. However, for the thermal conductiv- 
ity, k= n e + ki, we can only estimate the electronic 



part, which becomes more important in the itinerant 
electron limit. In the localized limit, k\ (the compo- 
nent due to the lattice) will be the only important 
contribution. Our results obtained from electronic 
structure calculations will be an upper limit for the 
overall zT but still helpful to understand how to 
optimize and engineer a better response in this and 
other related systems. 

The first parameter to analyze the TE response of 
the material is the Seebeck coefficient, that can be 
calculated independently of the scattering time in 
the constant scattering time approximation we are 
using for our calculations. Figure [3] shows the ther- 
mopower calculated for three values of hole doping in 
La2Ni04+5: 8= 0.05, 0.10 and 0.15, for several val- 
ues of the lattice parameter from 3.74 to 4.11 A, that 
simulate both the compressive and tensile strain lim- 
its (let us recall the unstrained value of the in-plane 
lattice parameter is 3.89 A). Results are presented 
using VCA (lower panel) for simulating doping or 
just displacing the chemical potential for the cal- 
culations in the undoped compound (upper panel). 
The values of doping chosen are, as explained above, 
in the range where the AF ordering survives and 
our calculations are more reliable, and also the in- 
terval where it was predicted the enhancement in 
TE figure of merit with doping for the unstrained 
compound^ 3 It was seen in previous studies 14 that 
large thermopower occurs for 5= and then gets 
rapidly reduced as doping is introduced. Here we 



5 




Temperature (K) Temperature (K) Temperature (K) 




Temperature (K) Temperature (K) Temperature (K) 



FIG. 3: (Color online.) Thermoelectric power of La2Ni04+^ as a function of temperature for three doping values: 5= 
0.05, (5= 0.10 and (5= 0.15 for various values of the in-plane lattice parameter (different strains). The upper panel was 
calculated for an undoped La2Ni04 shifting the chemical potential to the desired carrier concentration, whereas the 
lower panel was calculated using the VCA. The latter yields a more reliable and accurate description of the system. 
The exchange-correlation potential used was TB-mBJGGA. Observe that for various lattice parameters in the doping 
range 5= 0.05 - 0.10, thermopower exceeding 150 fiV/K are predicted. 



can analyze in more detail the effects of such doping 
levels in the thermopower together with the effect of 
strain. 

Looking at Fig. [3j one can see that in the ten- 
sile strain region (a > 3.90 A) the thermopower is 
enhanced, whereas in the compressive strain region 
it is reduced. This is the expected result because, 
as we have discussed above, enlarging the xy-plane 
leads to a wider gap and hence the thermopower will 
increase. We see that, as doping is increased, things 
get slightly more complicated than that, at the same 
time that thermopower is reduced. Putting every- 
thing together, values of the thermopower at room 
temperature in the range of S ~ 150 /iV/K are pre- 
dicted up to 5= 0.10 in the tensile strain regime, and 
up to 200 /xV/K for 5= 0.05. We observe in the lower 
panel of Fig. [3] that the results using the VCA dif- 
fer substantially in quantitative terms (some vague 
trends are maintained, but the values are very differ- 
ent to those obtained shifting the chemical potential 
in the undoped compound). We argue that using the 



VCA here, for these substantial doping levels, is nec- 
essary to yield a correct description of the system. 
We can see that the two methods give more simi- 
lar accounts at low-doping (S ~ 0.05), where shift- 
ing the chemical potential in the calculation for the 
undoped compound is enough to describe more or 
less accurately the system (it will work even better 
for S < 0.05), but going to larger doping levels, one 
needs to include the change in the overall potential 
caused by the oxygen excess in the system to get a 
good description. In the VCA-based calculations, we 
see that as doping increases, the values of the ther- 
mopower decay quite fast, except at high temper- 
atures. But even at high temperatures, the lattice 
parameter window where the response is promising 
becomes narrower, around 4.00 A. Focusing on that 
tensile strain value, we see that the values of the 
thermopower at room temperature are not promis- 
ing: 75 /iV/K at 5= 0.05, 40 /iV/K at 5= 0.10 and 
only 25 /iV/K at 5= 0.15, yet the high temperature 
response becomes enhanced. It has been observed 
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experimentally that La2Ni04+5 has an increasing 
thermopower at high temperatures 

However, we know that the larger the a parame- 
ter is, the wider the band gap gets (and that is why 
even upon doping large thermopower values can be 
obtained) but the conductivity will be reduced (in- 
plane hopping decreases). We can compile both ef- 
fects in the calculation of the dimensionless TE fig- 
ure of merit, that is presented in Fig. 3] for several 
doping levels and lattice parameters (the same ones 
we have used for the thermopower) and with the 
same calculation schemes described above. We see 
at low-doping, around 5= 0.05 and for a slight tensile 
strain (a > 3.90 A), values of the electronic-only TE 
figure of merit exceeding unity can be obtained even 
at room temperature for the larger lattice parame- 
ters tested, and at high temperatures for a very small 
tensile strain. When moving the system towards the 
compressive strain limit, even though in principle 
it should increase the conductivity, the gain through 
that term does not win over the loss in thermopower 
(which is substantial, as we saw above). We can 
see again that neglecting VCA and just shifting the 
chemical potential in the undoped compound will 
produce a large overestimation of the TE figure of 
merit (but the agreement with VCA is again bet- 
ter at small doping). As doping is increased, as we 
saw for the thermopower, a promising TE response 
only occurs at large temperature and for a very crit- 
ical range of in-plane lattice parameters (around a= 
4.00 A for 5= 0.10 and a = 4.11 A for 5= 0.15). 
This means that one has to optimize doping for each 
strain analyzed to obtain a significant TE figure of 
merit. For the larger tensile strain considered (4.11 
A), the behavior is quite peculiar, not being mono- 
tonic with doping. Similar non-monotonic behav- 
ior is observed also for the thermopower. As a in- 
creases, the x 2 —y 2 band gets lower in energy relative 
to the xz/yz bands, and eventually it gets so close 
to them that hole-doping introduces states of those 
bands very close to the Fermi level. The relevant 
contributions from those many bands make the be- 
havior and temperature evolution in that case much 
more complicated to analyze and understand. The 
large lattice parameter window with electronic-only 
figure of merit exceeding unity in the tensile strain 
limit only occurs for small values of doping around 
5= 0.05. 

As we have been arguing, one possible route to 
reduce the thermal conductivity is the growth in 



thin films. It has been recently observed how oxy- 
gen non-stoichiometries can lead to a significant re- 
duction of the thermal conductivity by creating new 
scattering centers that would hamper phonon prop- 
agation in SrTi03_5 thin films at very low doping 
levels i 4Q i 41 Experiments are required to understand 
if the same mechanism is at play in La2Ni04+5. If 
that is the case, the increase in the electronic fig- 
ure of merit via hole-doping could be a promising 
route to an enhanced TE performance in La2Ni04+5 
thin films. Moreover, the growth in the form of 
thin films will increase the electronic conductivity 
as has been shown in the past, 15 due to the larger 
weight of the in-plane conduction (the largely lay- 
ered structure impedes conductivity in the off-plane 
direction). Also, a thin film geometry is expected to 
reduce thermal conductivity by preventing the prop- 
agation of long- wave length phonons, apart from the 
possible reduction imposed by vacancy-related scat- 
tering. As we have discussed before, effects related 
to vacancy localization beyond the average treat- 
ment of the dopants we are considering here (both 
using VCA and a shift in the chemical potential) 
are beyond the scope of this work. The recipe, ac- 
cording to our calculations, would be the growth of 
La2Ni04+5 thin films with oxygen excess on the or- 
der of S ~ 0.10, grown on top of a substrate that pro- 
vides some tensile strain, and studied at high tem- 
peratures, where the most plausible applications of 
this oxide and the larger figure of merit will occur. 

To further validate our results, we can give an 
estimate of the overall zT (up to now we have 
studied the electronic part only). Considering the 
Wiedemann-Franz relation, the figure of merit can 
be rewritten as zT= <jTS 2 / n= n e S 2 / ' kLq, Lq being 
the Lorenz number with a value for free electrons 
L = 2.45 x 10" 8 (V/K) 2 and n < K e in general, 
but approaching 1 at large doping. At high enough 
temperatures, the lattice term typically decreases as 
1/T. This type of behavior can be observed in Ref. 
Il7l . Taking the thermal conductivity experimental 
values for the unstrained compound at the transition 
temperature from Ref. [I?], the electrical conductiv- 
ity from Ref. 15, and using our calculated values for 
the thermpower, we can obtain an estimation of the 
values of the overall zT. At 400 K, for a= 4.11 A, a 
zT of 0.3 can be achieved for 5= 0.05 while for 5= 
0.15 only a figure of merit of 0.04 can be obtained. 
For 5= 0.10 and a= 4.00 A, a lower value of 0.02 
is estimated at 400 K. However, at high temper a- 
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Temperature (K) Temperature (K) Temperature (K) 




Temperature (K) Temperature (K) Temperature (K) 

FIG. 4: (Color online.) Figure of merit of La2Ni04+<$ as a function of temperature for three doping values: 5= 
0.05, 5= 0.10 and 5= 0.15 for various values of the in-plane lattice parameter (different strains). The upper was 
calculated for an undoped La2Ni04 shifting the chemical potential to the desired carrier concentration, whereas the 
lower panel was calculated using the VCA. The latter yields a more reliable and accurate description of the system. 
The exchange-correlation potential used for calculating all the transport properties was TB-mBJGGA. Observe that 
for various lattice parameters in the doping range 5= 0.05 - 0.10, the electronic figure of merit exceeds unity. 



tures (above 700 K), zTs close to one (~ 0.8) can 
be obtained for a— 4.00 A and values of doping be- 
tween 0.05 and 0.10. These zT values are already 
very promising but if, in addition, the thermal con- 
ductivity of the system can be reduced with respect 
to the bulk value when grown in thin-film geometry, 
our calculations suggest that values of zT beyond 
one could in principle be achieved. 

IV. SUMMARY 

Our ab initio calculations for the compound 
indicate some promising features in hole-doped 
La2Ni04 + 5 as a possible oxide thermoelectric ma- 
terial. If the thermal conductivity can be reduced 
by nanostructuring, e.g. in the form of thin films, 
the system could show enhanced thermoelectric per- 
formance at low hole-doping levels, attainable by 
the appropriate control of the apical oxygen con- 
tent. Our calculations show that a region with rel- 



atively large Seebeck coefficient exists in this com- 
pound at small doping levels and the appropriate 
tensile strain, within the realistic AF description. A 
careful experimental study needs to be performed 
in this respect controlling the oxygen content, and 
also making thin films with the appropriate oxygen 
composition. Hole-doping will increase the electrical 
conductivity, as the thin film geometry does, which 
together with the reduction in thermal conductivity 
could leave room for an improvement of the thermo- 
electric figure of merit of this layered nickelate. We 
suggest to explore the region of tensile strain, grow- 
ing on substrates with a lattice parameter slightly 
larger than that of SrTi0 3 (about 3.95 - 4.05 A). 
Our calculations predict a reasonable thermoelectic 
performance (with figure of merit close to 1, zT ~ 
0.8) in those conditions of doping and tensile strain 
at high temperatures. 
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